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Abstract 
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Introduction 



New interesting solutions for gravitational instantons exist in the literature. One of them is the helieoid 
solution of Nutku [1,2 ]. Since this solution has eurvature singularities, it has not been studied extensively 
aside from three articles [3,4,5 ]. These three papers study the solutions of the Dirac equation using this 
solution as a background. Our work on these solutions can be extended by studying them in a manifold with 
a boundary. Taking a domain that excludes the origin and infinity for the radial variable cures the problems 
associated with the infinities of these solutions. 

One of us (M.H.), with collaborators, had studied the related problems in the presence of an instanton and a 
meron in the past [6,7 ]. Similar problems have been studied recently for a spherical bag [8 ]. In these three 
papers, to conserve self adjointness, as well as chirality and charge conjugation, spectral boundary conditions 
[9,10 ] were used. The spectral boundary conditions were also extended to the Dirac equation with torsion 
[11]. 

Here we study the solutions to the Dirac equation in the background of the Nutku helieoid metric. We 
use Euclidean signature for our metric. This makes our differential operator of the elliptic type. We find 
that although some of the solutions diverge at the origin, they are normalizable when integrated with the 
measure , where r is the volume. We also see that the solution diverges as x goes to infinity, since they 

are expressed in terms of modified Mathieu functions. Forcing them to be regular at the origin makes them 
divergent at infinity. Furthermore, the metric used has curvature singularities at the origin. Although one 
can define the index on certain noncompact manifolds, the conditions of "bounded geometry" are needed 
[12,13 ]. These are the conditions we do not meet. Thus we have to restrict our space at a boundary. This 
necessitates specifying the boundary conditions. 

We know that local boundary conditions can be used in even dimensions, although, it was shown [6 ] that 
the spectral boundary conditions [9,10 ] are the only self adjoint one that also conserves the 7^ and the 
charge conjugation symmetry of the Dirac operator. The chiral Dirac operator, however, requires spectral 
boundary conditions in even dimensional spaces [14]. The non-local, spectral boundary conditions should 
be used for the odd dimensional case due to obstructions that were first pointed out by Atiyah and Bott 
[15 ]. Here we will study the problem both in odd and even dimensions. We will use the spectral boundary 
conditions in both cases. 

The non-local spectral boundary conditions were introduced by Atiyah- Patodi-Singer while they were inves- 
tigating the Hirzebruch signature theorem for manifolds with boundaries [15,16,17]. When a manifold has a 
boundary, we may use local elliptic boundary conditions, like those of Dirichlet, Neumann or Robin, if the 
manifold has even dimensions. We may be sacrificing the symmetries cited above, though. If the manifold is 
odd dimensional, however, we have to use the non-local, spectral boundary conditions. Atiyah and Bott [15 ] 
found out that in general there are topological obstructions to finding acceptable local boundary conditions. 
As described in [15,16 ] the crucial point is defining an acceptable elliptic boundary problem for the signature 
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operator whose index is the signature of the manifold. We call a complex elliptic if we can define an elliptic 
operator on it. The exterior algebra can be split into two distinct elliptic complexes. The first is the de Rham 
complex, which is related to the Euler characteristic. The second is the signature complex [17 ]. Although 
the de Rham complex admits local boundary conditions, it can be shown that for the signature complex, as 
well as for the spin and the Dolbeaux complexes, there does not exist a local boundary condition so that this 
complex with the local boundary condition is elliptic [18 ]. In the first order case, like the Dirac operator, 
there is a natural pseudo-differential operator given by the projection on the generalized eigenvectors with 
eigenvalues with plus or minus real parts. This operator leads to a well posed boundary value problem for 
the signature complex ( as well as the spin and the Dolbeaux complexes.) This fact requires the use of the 
spectral boundary conditions, which are non-local. 

In defining the index of the operator, the presence of a boundary necessitates extra terms in addition to the 
Euler number. For the case without a boundary, knowing only the curvature suffices to compute the index 
[17 ]. One calculates the integral over the whole manifold of the same characteristic classes. If the manifold 
has a boundary, then, we have to add a term obtained from the integration of the Chern-Simons form, which 
is written in terms of the connection, the curvature and the second fundamental form determined from the 
normal to the boundary. Another term that should be added is the r] invariant, which is calculated from the 
eigenvalues of the operator restricted to the boundary. We also add a term proportional to the number of 
zero eigenvalue solutions of this operator. The index is equal to the sum of these four terms. 

Here wc study the Dirac equation in the background of the hclicoid metric in five and four dimensions in 
consecutive sections. We do not attempt to calculate the index for our Dirac operator. Our only aim is to 
make a sense out of this incomplete metric by restricting the domain and using the appropriate boundary 
conditions for this problem. 



Solutions of Dirac Equations 

The relevance of the helicoid- catenoid metric, as a simplest minimal surface in Euclidean space resulting 
in a gravitational instanton [1], for certain cosmological models is given in a recent preprint [21]. One can 
transform from the helicoid to the catenoid metric by a simple transformation [2]. The importance of five 
dimensional spaces is also stressed in the same paper. The helicoid metric is also mentioned among the 
metrics related to M theory by Gibbons et al [22]. This metric is also studied by Valent in [23,24,25]. These 
are other reasons why we try to give a meaning to the solutions of Dirac equation in the background of the 
Nutku helicoid solution both in four and five dimensions. Here we briefly introduce the metric and write 
formally the Dirac equation in the background of this metric. 

The Nutku helicoid metric is written for the four dimensional Euclidean space. It is given as [2] 
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ds^ = ^ " J dr^ + (r^ + o2)d02 + ( 1 + ^sin'^e) dy^ - ^sin20dydz+ (l + ^cos^9 ] dz% [1] 



where 0<r<oo, 0<^< 27r, y and z are along the KiUing directions and will be taken to be periodic 
coordinates on a 2-torus. This is an example of a multi-center metric. This metric reduces to the flat metric 
if we take a = 0. 

If we make the following transformation r = asinhx, the metric is written as 

a2 2 
dsA = —sinh2x(dx^ + dO^) + — — -—Usinh?x + sin^9)dy^ — sin20dydz + (sinh^x + cos^6)dz^]. [21 
2 sinh2x 

We use the NP formalism [26,27] in four Euclidean dimensions [28,29,30 ]. The details of how the tetrad is 
chosen, how the 7** matrices are formed for this explicit case are given in [5 ]. 
We use the Dirac operator i'y^V^ where 

= - r^. [3] 

The 7 matrices satisfy 

{r,Y} = 2g^\ [4] 
The spin connection is written as = j7?]Ii7i/. 

We can trivially extend this structure to five dimensions The addition of the Euclidean "time component" 
to the previous metric gives: 

rfs^ = dt^ + dsl. [5] 

Going to five dimensions does not increase the number of equations we have to solve, since in both four and 
five dimensions we can use four component Dirac spinors. We can not interpret the fifth dimension as time, 
since this will bring problems with causality. If we use non-local boundary conditions, these conditions are 
set for all times. Abrikosov has devised a new method [31 ] which cures the causality problem. We could not, 
however, generalize this method to our case, since our little Dirac equation still couples three components. 
We, therefore, in both even and odd dimensional cases, study the solution with the Euclidean metric. We 
include both cases to compare and contrast the similarities and difii'erences of these closely related cases. 

Solutions in five dimensions 

We start by studying the solutions to the Dirac equation in the background of the Nutku helicoid metric in 
five dimensions. Our motivation is to see whether one can find a domain where solutions in the background 
of this singular metric can be defined. We will try to show that one can obtain a satisfactory solution to this 
problem. 



We write the system in the form Lip = Aip and try to obtain the solutions for the different components. 
Our aim is to write the upper components in terms of the lower components. We will impose the boundary 
conditions on the upper components in terms of derivatives of the lower components below. 
The equations read 



■\/2 dv sifih^x 

-.{{d^+ide)'^z +a[cos{e + ix)dy + sin{9 + ix)d^\^4, —= ^4*1} = A*i, [6a] 



^ -.{{d^-ide)^4 -a[cos{e-ix)dy + sin{e-ix)d,]^3-^^^^^^^^^dt^2} = A^2, [66] 



aV sinh2x \/2 

-.{{d^-id0+coth2x)^i - a[cos{9 + ix)dy + sin{e + ix)d:,]'^2 + ^ 9**3} = A*3, [6c] 



aV sinh2x V2 

-.{{d^ + ide + coth2x)'^2 + a[cos{9 - ix)dy + sin{e - ix)d:,]'^i -\ p ^4*4} = A*4. [Qd\ 



aV sinh2x ^/2 
These arc coupled equations involving three different components. The usual method to obtain solutions at 
this stage is the separation of variables method. Three of the variables, t, y, z, define the Killing directions. 
This is exploited by writing the solution as a product of exponentials in y, z, t times a function of x and 0. 

We take ky = kcos{^), kz = ksin{(p), since then we can absorb the variable <p in functions of the remaining 
variables x and 0. 

We note that these four equations are not similar in form. The transformation ^12= / ^ , fi 2 is used 

' V S2nh2x ' 

for the upper components. This transformation eliminates the coth2x terms in the third and the fourth 
equations. With this transformation, the equations read 

{idx+ide)'S3 +iak[cos{9- (t> + ix)]'i'i-iaktfi} = A-^fi, [8a] 

V2 



{id^-id0)^4 -iak[ms{9-(f)-ix)]^3-iaktf2} = A-^f2, [86] 

v2 

{-d^ + ide)h + iak[cos{9 - + ix)]f2 + iakt^^ = -A ^^'^/^^^ ^s, [8c] 

{-dx - ide)f2 - iak[cos{9 - ^ - ix)]fi + iakt^i = -A — ^-4. [Sd] 

\ z 

We solve our equations in terms of /i,2 and substitute these expressions in equations, given above. This 
substitution gives us second order, but uncoupled equations for the lower components. 
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{dxx + dee + ^[A;^ (-cos[2(6' - (j))] - cosh2x) - {kl + A^)sinh2x\^ *3,4 = 0. 



[9] 



We can separate this equation into two ordinary differential equations by the ansatz ^'3^4 = R{x)S{6 — (j)). 
Separation of equation(9) gives us two ordinary differential equations. The equation for S reads 

deQS{Q) - fc2cos(2e) - S{Q) = 0, [10] 

where {6 — (p) = 6. This equation is of the Mathieu type and the solution can be written immediately. 

S{e) = CiSe{n, -<!>)+ C2So{n, 6 - c/,) [11] 

The solutions should bo periodic in the angular variable O. This fact forces n, the separation constant, 
to take discrete vahies. It is known that the angular Mathieu functions satisfy an orthogonality relation 
such that functions with different n values are perpendicular to each other. Here, we integrate the angular 
variable from zero to 2tt. One can fix the normalization constant according to the chosen normalization, 
whether it is according to the McLachlan or the Morse-Stratton convention [32 ]. 
The equation for R{x) reads 



d^x - [Yik'^cosh2x + {kl + K^)sinh2x] + n \ R{x) = 0. [12] 
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whose solution can be reduced to the form 

R{x) = DiSe{n, Ag, i{x + b)) + D2So{n, Aq, i{x + b)). [13] 

Here Ci, C2, fi, £>2 are arbitrary constants. The other constants will be defined below as we explain how 
one can reduce our initial equation to give these solutions. The solutions for the lower components il>3,tl>4 
are given in terms of sums over n and integrals over k,(p,kt we used in the separation ansatz. We think it 
is amusing to explain how these solutions are obtained. To get this solutions one has to go through several 
steps. 

In reference books the solution of eq. (13) is not listed as a Mathieu function. In fact, we find that the 
solution of this equation is expressed in terms of double confluent Heun functions [33 ]. We denote these 
functions as Hn in our expressions. 

R{x) =CiHdI 0, — n, a^(A^ + k^), + n, tanhx j 

+C2HD ( 0, — n, a^{h? + k^), + n, tanhx \ 

[14] 



/—dx 
Hd (0, ^ - n, a2(A2 + fc^), ^ 



+ n, tanhx) ^ 



Normally one takes the first function and discards the second solution. 

Reading through literature [34 ] we suspect that the solution we found can be expressed in terms of Mathieu 
functions after performing proper transformations. We can show that we can indeed express the result in 



terms of Mathieu functions by using a simpler independent variable than the one given in eq. (14) after 
performing few transformations. We define 

A2 = -n [16] 



A3 = - 

and use the transformation 



^3 = -^, [17] 



z = [18] 
Then the differential operator in eq. (12) is expressed as 

4z'd,, + 4.5, + {^^^)z + A, + (^^4^)^) / = 0. [19] 

This equation is still of the double confluent Heun form, since it still has two irregular singularities at zero 
and infinity. We define A4, = (^^^^) and A5 = i^^^). If we take 

J^u = z, [20] 
V ^5 



and 



and set Aq = ^/A^As we get, 



-=i(«+^) [21] 



(u;2 - l)d^^ + wdy, + i^w + / = 0. [22] 

With the new transformations, we have traded the irregular singularity at zero by two regular singularities 
at plus and minus one. This is the same singularity structure of the Mathieu equation. The solution of this 
equation is indeed expressible in terms of Mathieu functions. It is given as: 



R{'w) = Se l^n, Aq, arccosy "~"^~~"J + '^'o yn, Ae, arccos\j ' [^3] 

where A(, = ^[/c^ — (fcj + A^)^]^/^. Going back through the transformations we made, it is not hard to 
express arccos'^ ^£±1 in terms of our original variable a; up to a constant, i{x-\-h). We simply write 



-2x 



I fc2_fc2_A2 



cosh'^{x + h). [24] 



Here 



e 



-26 _ / fc^ kj -^^ rr,rl 
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Taking the arccos of this expression gives the result given in eq. (13). 

After these transformations, we see that the solution of the lower components of the Dirac equation can be 
expressed in terms of functions that are regular at zero. We can not say this for the upper components, 
though. They are expressed in terms of these solutions and their derivatives divided by a function, which 
blows up at zero. Even if we take the so called odd Mathieu functions, which can be expanded in terms 
of hyperbolic sine functions, their derivatives will be hyperbolic cosine functions. To obtain the upper 



components we have to divide them by V sinh2x . Then these functions will still blow at the origin. 
A finite scalar product can be defined around the origin for these solutions in the form 

j %-^i^dT,i=l-A, [26] 

in a finite domain, including the origin. Here dr is the volume element in our five dimensional space. 
Repeated indices are not summed over. ^ is the square root of the determinant of the metric, necessary 
to get an invariant volume element. The zero of the invariant measure cancels the singularity of the wave 
functions at the origin. 

We find, however, that the solutions are not normalizable as x goes to infinity. Furthermore, our metric has 
curvature singularities at the origin. Although, using our new measure, we can make all four components 
normalizable at the origin, we still define our solutions for the domain < a; < F. We integrate x variable 
in the domain up to the point x < F, where F is the point the function starts to diverge. Since our radial 
solution is multiplied by the angular solution and the exponential function to make up the total solution of 
the Dirac equation, the orthogonality of the angular solutions for different values of the discrete n makes our 
solutions orthogonal to each other. By dividing by the appropriate factors, we can normalize them. 
We note that the domain where the solutions are normalizable is restricted. These differential equations do 
not have a meaning unless we define the boundary conditions the solutions obey at the boundary of our 
domain. Since our system exists in an odd dimensional manifold with a boundary, we have to study this 
equation using spectral boundary conditions [9,10 ] . 

The method used in applying these boundary conditions requires first studying the little Dirac equation, the 
tangential operator of the Dirac operator restricted to the boundary, where the variable x takes a fixed value 
xo- We have to take xo greater than zero since the second fundamental form defined by our choice of the 
normal to the boundary diverges at a;o = [18 ]. 

For this purpose we write the equations in the form given in equations (8), i.e. after ^'1^2 are transformed 

to /l,2. 

\/2 ink* 

— {ide^3 + iakcos{e - + ixo)^i - -^h} = Xfi, [27a] 
a y 2 

\/2 inkt 

—{-ide^i - iakcos{6 - <^ - ixo)^^ - -^h} = A/2, [276] 
a y 2 

\/2 inkt 

— {-ideh -iakcos{9-(l> + ixo)f2 + ^'^3} = X^3, [27c] 
a y 2 



\/2 inkt 

— {idef2 +iakcos{e-(f>-ixo)fi + ^^4} = \^4- [27d] 
Here A is the eigenvalue of the little Dirac equation. 

We could not obtain analytical solutions of those equations in terms of known functions. Wo could not even 
write uncoupled equations in the sc;c;ond order. One needs to go to fourth order in derivatives to be able to 
write equations that involves a single unknown function. 

At this point wo follow closely our references [6,7,8,31]. We formally expand our solutions at the boundary, 
fixed by two values of xq in terms of eigenfunctions of the little Dirac equations with both positive and 
negative eigenvalues A. For the lower components we have 

*3,4(0,aro) = ^/iA(e,aro) [28] 

for fixed values of kt,ky,kz. We set 

*3,4(0,a;o)|oB = Y,hx{Q,xo). [29] 

A>0 

The negative A eigenvectors are all set to be zero at the boundary. 

Then, we solve /i^2 in terms of ^'3.4 using the equations (8a, 8b) and fix the x values to xq on the boundary. 
Note that since x = .xo, the derivative with respect to x is evaluated at this point as well as the terms 
without the x derivative. We can, in general, use the expansion given for "^3^4 on the boundary, in terms of 
its eigenfunctions. 

*3,4(0,a;o) =^/iA,3,4(e,a;o). [30] 

A 

This sum is over all values of A. In fixing the values of /i,2 in terms of ^'3,4 on the boundary, we use only 
part of the expansion of ^'3 4(6, x) where 

*^,4(e, x)\aB = Y^ hxM^o, e). [31] 

A<0 

In other words, we write ^'^2(01 ^)I<3B using the expressions obtained from /i^2 in terms of the negative A 
values of 'i'^ ^{Q, x)\dB ■ These boundary conditions arc non local, but are shown to be the only consistent 
ones for odd dimensional Euclidean spaces by Atiyah-Patodi-Singer. 



Solutions in four dimensions 

Here we repeat the calculations given above after setting kt to zero. Our aim is to show that the five 
dimensional case, studied above, is not essentially diff'erent from the four dimensional case. Since both 
solutions can be expressed in terms of Mathieu functions we can compare the four and five dimensional 
cases. These equations with A = were studied in ref. [5 ] . The solutions obtained in this reference for the 
four dimensional case can be summed to give the Green's function similar to the calculation done in ref. [2 
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]. This is not possible for the five dimensional solution of ref. [5 ]. For the solutions with A ^ both in 
five and four dimensions this property does not exist. The radial and the angular equations have different 
constants which does not allow to obtain the Green's function of the equation by using summation formulae 
of Mathieu functions [35]. 

Our solutions exist only in a finite domain for the variable x, hence we have to use the appropriate boundary 
conditions at this point. Using the orthogonality of the angular Mathieu functions for different values of the 
discrete parameter n, we can show the orthogonality of these solutions. 

At this point, although we can use local boundary conditions in this case, we choose to use the spectral 
boundary conditions of Atiyah-Patodi-Singer to conserve chirality and charge conjugation. We are keen not 
to break chirality by hand, since in the standard model this symmetry is expected to break spontaneously, 
resulting in confinement. There are cases, however, to describe the correct physics, where one may want to 
break this symmetry to force the system to one of the broken phases. Using chirality breaking boundary 
conditions may be one way to achieve this task. This is discussed in detail in reference [36]. In this work, 
we do not want to break the chiral symmetry of the model. We, therefore, adhere to the spectral boundary 
conditions. 

We write the system in the form = where is a four component spinor, and try to obtain the 
solutions for the different components. Then our equations are similar to the ones given in the equations 
(6). The only difference is taking kf equal to zero. We see that three components are still coupled in our 
equations. 

The method of solution is exactly like the one used in the previous section. To get our solutions we use the 
separation of variables method. We write the solution as a product of exponentials in y, z times a function 
of X and Q. 

^. = ei(fc„y+fc.^)^.(a.^^)_ [32] 

The same transformations are used as those in the five dimensional case to reduce 5'i_2 to /i^2- We solve 
our equations in terms of /i^2 and substitute these expressions in equations, given above. We end up with 
second order, but uncoupled equations for the lower components. 

[dxx + See - Y [fc^ {cos2{Q - 0) + coshflx) + IS?sinh2x\^ ^3,4 = 0. [33] 

We can separate this equation into two ordinary differential equations by the ansatz ^'3^4 = R{x)S{6). For 
S{Q) we get an equation of the Mathieu type and the solution can be written immediately. 

S{e) =CiSe[n,^,0-(P]+C2So{n,^,e-(P). [34] 

Solution for R{x) can be reduced to 

R{x) = DiSo{n, B, i{x + h')) + D2Se{n, B, i{x + b')), [35] 
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Here Ci, C2, -Di, -D2 are arbitrary constants. B = . b' is defined as in eq.(25) with kt = Just note 

that n is the separation constant which has to take discrete values to get a periodic solutions for the angular 
Mathieu equation <S'(0). The solutions for the lower components ^'3,^4 are given in terms of sums over n 
and integrals over k,(f) we used in the separation ansatz. 

We find that we have the upper solutions, 1,2 are divergent at the origin, whereas the lower ones are finite. 
Both of our solutions diverge at infinity. For the same reasons as given for the five dimensional case, we have 
to limit the domain of our solutions at two finite values of x. 

To impose these boundary conditions we need to write the little Dirac equation, the Dirac equation restricted 
to the boundary, where the variable x takes a fixed value xq. We choose to write the equations in the form, 

/2 

— {i5e*3 +ikacos{0 - (j) + ixo)'^4} = >^fi, [36a] 

^{-id0^4 - iakcos{e - ^ - ixo)^^] = A/2, [366] 
\/2 

— {{-ide)fi - iakcos{0 - + ixo)f2} = A*3, [36c] 
v/2 

-^{{ide)f2 +iakcos{e-(l)-ixo)h} = \^i. [36d] 
Here A is the eigenvalue of the little Dirac equation. 

We could not obtain analytical solutions of these equations in terms of known functions. This is the same 
result as in the five dimensional case. For us it seems very curious being able to solve similar system 
of partial differential equations, but not even being able to decouple them when this system reduces to 
ordinary differential equations on the boundary. One possible explanation is that 6 — (f>±ix act as z and z of 
complex variables. Sometimes it is easier to find functions of this pair as solutions is easier than a function 
of a single real variable. Actually for the full Dirac equation, Sucu and Unal [3] find solutions in a closed 
form. Same technique, however, does not seem to work when x = xq. 

From equations (36) we see that we can take ^'4 as the complex conjugate of and /2 as the complex 
conjugate of /i. We take /i ~ pi + iqi and 4*3 = ps + iqs- We end up with coupled differential equations for 
these functions. One notes, however, that separating the real and imaginary parts of each solution, we get 

deePi + ak[—sinQcoshxo — {sin^Qsinh^xo + cos^Qcosh^xo)]pi — akcosQsinhxoqi = X^pi, [37] 



deeQi + ak[sinQcoshxo — {sin^Qsinh^xo + cos^idcosh^xo)]qi — akcosQsinhxopi = X^qi. [38] 

We see that the eigenvalue of the little Dirac equation comes only quadratically, showing a symmetry for its 
positive and negative values. This will make the r] invariant, defined as 

r] = HTOs^oy^szgn(Ai)|Ai|~^, [39] 
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which is needed for an index calculation, zero. We also see a symmetry between pi and qi , namely pi (©) = 
©)• If we, instead, eliminate pi,qi and write our equations for pa, 93, we get exactly the same equations. 
We expand our solutions at the boundary, fixed by two values of xq in terms of eigenfunctions of the little 
Dirac equations with both positive and negative eigenvalues A. 

*f(©,a;o)=^9i,A(e,a;o) [40] 

A 

for fixed values of ky,kz. We set 

*3,4(0, ^)\aB = ^ <?A,3,4(e, Xo). [41] 
A>0 

The negative A eigenvectors are all set to be zero at the boundary. 

The boundary conditions on the upper components are imposed exactly in the same manner as explained 
in the five dimensional case, namely we solve for /i,2 in terms of ^'3,4 using the equations 

^{{d^+idg)^3 +a[cos{e-cl> + ix)]^i} = Afu [42] 
v/2 

^{{d^-ide)'^4 -a[cos{e-(j>-ix)]^3} = Af2, [43] 
and fix the x values to xq on the boundary. 

We can in general use the expansion given for ^'3^4 on the boundary, in terms of its eigenfunctions, eqn. 
(40). This sum is over all values of A. In fixing the values of /i,2 in terms of ^'3,4 on the boundary, we use 
only the part where A < 0. These boundary conditions are non local, but they respect self adjointness and 
conserve 7^ and charge conjugation symmetry. 

Conclusion 

Here we tried to give solutions of the Dirac equation in five dimensions for the Nutku helicoid metric in a 
bounded region. We found out that they can be reduced to Mathieu functions, which is also the case in four 
dimensions. We imposed formally [8,31] the non-local spectral boundary conditions of Atiyah-Patodi and 
Singer [9,10,6,7] on these solutions, which are the only correct boundary conditions in odd dimensions. Our 
main goal in this paper is to define the solutions of the Dirac equation in the background of a incomplete 
metric, namely the Nutku helicoid solution consistently. 

A related work would be to calculate the index of the Dirac operator in four dimensions in this background. 
Prom our solutions we can calculate both the bulk and the surface term easily. From the form of the little 
Dirac equation, we see that the 77 invariant is zero. We could not obtain analytical solutions of the little Dirac 
equation, though. We give examples of the numerical solutions we found of the zero eigenvalue equations of 
the little Dirac equation in Figure 1. These pictures do not correspond to any of the functions we encounter 
in the literature [32,33,35,37]. Therefore, we know that the zero mode solutions of the little Dirac equation 
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exist. We, however, do not know their analytical expression, hence their number. As a result, we could not 
calculate the index in this paper. 
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FIGURE CAPTIONS 

Figure 1 : Numerical solution for zero mode solutions of pi and qi in four dimensions for xq = 0.005, a = 
l,k= 1. 
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Figure 1 
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